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ABSTRACT 

The collapse of weakly turbulent prestellar cores is a critical stage in the process 
of star formation. Being highly non-linear and stochastic, the outcome of collapse 
can only be explored theoretically by performing large ensembles of numerical 
simulations. Standard practice is to quantify the initial turbulent velocity field in 
a core in terms of the amount of turbulent energy (or some equivalent) and the 
exponent in the power spectrum (n = — <i log P^/d log fc) . In this paper, we present a 
numerical study of the influence of the details of the turbulent velocity field on the 
collapse of an isolated, weakly turbulent, low-mass prestellar core. We show that, 
as long as n > 3 (as is usually assumed), a more critical parameter than n is the 
maximum wavelength in the turbulent velocity field, A MAX . This is because A MAX 
carries most of the turbulent energy, and thereby influences both the amount and the 
spatial coherence of the angular momentum in the core. We show that the formation 
of dense filaments during collapse depends critically on A MAX , and we explain this 
finding using a force balance analysis. We also show that the core only has a high 
probability of fragmenting if A MAX > i? CORE /2 (where i? CORE is the core radius); 
that the dominant mode of fragmentation involves the formation and break-up of 
filaments; and that, although small protostellar disks (with radius R DISK < 20 AU) 
form routinely, more extended disks are rare. In turbulent, low-mass cores of the type 
we simulate here, the formation of large, fragmenting protostellar disks is suppressed 
by early fragmentation in the filaments. 
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1 INTRODUCTION 

Although there is strong observational evidence that circum- 
stellar disks are formed during the early (Class and Class 
I) phases of protostellar evolution, estimates of their masses 
and extents are uncertain by up to a factor of 10 ( J0r gensen| 
et al. 2009). Since protostellar disks that are sufficiently mas- 



sive and extended to fragment might be an important site for 
forming low-mass stars, brown dwarfs and planetary mass 
objects (e.g. Stamatellos et al.|2011 1, it is important to un- 
derstand the circumstances under which such disks can form. 

In this paper, we study the formation of protostellar 
disks in collapsing, weakly turbulent cores, and evaluate the 
influence of the scale of turbulent perturbations and the 
net core angular momentum. Q In contrast to models based 
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1 Note that we are here using the term turbulence in the loose 



on rigid rotation (|Walch et aL]|2009| |Machida et al.||2010" 
|Machida fc Matsumoto||2011[ ), it has never been demon- 
strated unambiguously that the net angular momentum in 



a turbulent core (Dib et al. 20101 significantly affects the 



size of the protostellar disk that it spawns. Rather, previous 



numerical models (e.g. Goodwin et al. 2004a Walch et al 



20101 find no correlation between core angular momentum 



and disk size, suggesting that disk size is determined by 



sense of random (and here statistically isotropic) macroscopic mo- 
tions over a range of length-scales, and not in the more restrictive 
sense of fully developed turbulence, in which energy is injected on 
large length-scales and cascades through a large inertial range 
of length-scales, before being dissipated on much smaller length- 
scales (e.g. |Kritsuk et al.|2007| . Here, turbulence is simply a device 
for seeding a star-forming cloud or core with the fluctuations that 
will eventually lead to fragmentation; the term is routinely used 
in this loose sense by those who simulate collapse and fragmen- 
tation. 



2 S. Walch, A. Whitworth, P. Girichidis 



details of the specific turbulent velocity field. However, pre- 
vious work dKlessen et al. 120001 |Fisher|2004 



Goodwin et al. 



2004a b Matz ner fc Levin||2005| IKratter fc Matzner 
Goodwin et al.||2006 



2006 



IKrumholz et al.| |2007| jAttwood et al. 



2009| |Walch et al.||2010| > has not investigated fully the pa- 
rameter space used to initialize turbulent prestellar cores. 

A random, Gaussian, turbulent velocity field is charac- 
terized - in a statistical sense - by five parameters, (i) Some 
measure of the total amount of turbulent energy, for exam- 
ple the root-mean-square turbulent velocity, or the 



mean turbulent Mach number, M. Observations (Goodman 



et al. 1998 Barranco & Goodman 1998) suggest that tur- 



bulent velocities in low-mass cores are approximately sonic 
(A4 ~ 1). (ii) The partition of energy between solenoidal 
and compressive modes (e.g. Federrath et al. 20081, where 



the statistical equilibrium distribution is 2:1 solenoidal to 
compressive, (iii) The slope of the turbulent velocity power 
spectrum n, where Pk oc k~ n , and n is typically between 
3 and 4 (Burkert & Bodenheimer 2000 note that here n 
is defined so that - for example - the Kolmogorov scal- 
ing index corresponds to n= 11/3). (iv) The wavelength of 
the largest turbulent perturbation, A MAX = 2_R CORE /fc MIN . 
(v) The wavelength of the smallest turbulent perturbation, 



A 



1R r 



Previous studies of low-mass core 



MIN 

collapse - and of star cluster formation from the collapse of 
larger, more turbulent molecular cloud cores, e.g. |Bate et al~1 
( 2003[ ); |Bonnell et al.| ( |2003[ >; |Batel ( |2009a|b| ) - have not ex- 
plicitly specified the last two parameters (A MAX ,A MIN ). In 
this paper we evaluate their influence in more detail. 

As long as (a) the power spectrum is sufficiently steep, 
n > 3, and (b) the inertial range of the initial turbulent ve- 
locity field, A MAX /A MIN = fc MA x/ fc MiN> is sufficiently large, 
then A MIN is unimportant, since very little turbulent power 
is invested in the shortest wavelengths and it is divided be- 
tween many modes. However, A MAX has a major impact on 
core collapse, fragmentation, and disk formation. For the low 
levels of turbulence typical of low-mass cores, dynamical fila- 
ment fragmentation requires A MAX > R CORE ; fragmentation 
is very rare when A MAX < J? CORE /2. In addition, core angu- 
lar momenta and the radii of protostellar disks both increase 
with increasing A MAX . These differences arise because large- 
scale turbulence promotes the formation of large, coherent 
filaments. Such filaments not only fragment, but also deliver 
streams of material with disparate specific angular momenta 
into the center of the core, where this material then forms 
large disks. 

The structure of the paper is as follows. In section 2 
we describe the numerical method. We present the results 
of our simulations in Section 3, discuss the results in section 
4, compare with previous work in section 5, and summarize 
the main conclusions in section 6. 



2 INITIAL CONDITIONS & NUMERICAL 
METHOD 

2.1 Turbulence 



these velocities result in a net angular momentum (Gold- 



Turbulent cores have irregular internal velocity fields (Bel- 
loche et al.|2001||Andre et al.|2007||Maruta et al.|2010[), and 



|smith fc Arquilla|1985l|Dubinski et al.|1995"l|Goodman et al 



1993 Barranco & Goodman 1998 Caselli et al. 20021. Ji 



jina et al. ( 1999 1 estimate that a low-mass core typically 
has a ratio of turbulent to gravitational energy in the range 
< 7turb ~ 0.5, and a mean specific angular mome ntum 



Burkert & Bodenheimer (20001 have 



shown that these features can be reproduced if the turbu- 
lence has a power spectrum of the form Pk oc k~ n with n = 3 
or n — 4, and this finding has been employed by several au- 



thors (e.g. Fisher 2004| 


Matzner & Levin||2005| Kratter & 


Matzner||2006| Krumho 


z et al.||2007| |Kratter et al.||2008| 


Walch et al.||2010 


1. 



Here, we create random, turbulent velocity fields with 
the same ansatz. First, we generate random Gaussian ve- 
locity fields with Pk oc fc -4 (Burgers turbulence), populat- 
ing the wavenumbers fc MIN ^ k ^ fc MAX in Fourier space. 
Next, we map the velocity fields onto a uniform 128 3 grid, 
and scale them to the required root-mean-square velocity 
(or equivalently, to the required Mach number, assuming a 
sound speed of 0.2 km s" 1 ). Finally we compute the initial 
velocities of individual SPH particles by linear interpolation 
on this grid. This setup results in a ratio of solenoidal to 
compressive modes close to 2:1. 

Previous simulations of turbulent cores have not re- 
ported the range of wave-numbers populated in the ini- 
tial turbulent velocity field. We demonstrate here that this 
is a fundamental issue. Specifically, the choice of fc MIN 
is critical, but the choice of fc MAX is unimportant, pro- 
vided fc MAX 3> k MIN . To show this, we consider fc MIN 



I ' :' 1, 2, 4, corresponding to A^ 



2R r 



4 - R CORE. S-RcORE. -RcORE. ^CORe/ 2 ' FOT eaCn V£llue ° f 

fc MIN we simulate five different realizations (by using five 
different random number seeds 100, 200, 300, 400 and 500), 
in order to obtain better statistics. In each case we adopt 
fc M4s = fc MIN + 8. Since n = 4, there is little power on the 
shorter wavelengths, and so the choice of fe MAX is unimpor- 
tant. 



2.2 Initial core properties 

We base the initial structure of the simulated core on the 
observed properties of the A-MM8 core in Ophiuchus. 



Using 850 (im observations, Simpson (20101 estimates 
the mass to be M core = 1.28 M , and the azimuthally av- 
eraged FWHM diameter to be D FWHM = 2000 AU (A- 
MM8 is slightly elongated, with aspect ratio < 1.2, but 
we neglect this). The core is modeled with the density pro- 
file of a critical Bonnor-Ebert sphere (i.e. the dimension- 
less boundary radius is £ B = 6.451), but it is too cold to 
be in hydrostatic equilibrium. The outer radius of the core 
is given by i? C ore = 1-25 Dfwhm = 2500 AU (on the as- 
sumption that the dust temperature is uniform and the 
850 /jm emission is optically thin). The central density is 



6.5 x 10" 
i.-,.- - 8.3 kyr. 



Pc 
t 



gem ; and the central free-fall time is 



Detailed 3D radiative transfer modeling of A-MM8 by 
Stamatellos et al. (2007) indicates that the mean tempera- 
ture is T — 11 K, and therefore the ratio of the thermal to 



gravitational energy is a, 



= 0.017. 



Andre et ah] ([2007] report a FWHM velocity width for 
the N 2 H (1 - 0) line of Au FWHM = 0.384kms _1 . If the 
non-thermal velocity dispersion is attributed to turbulence, 
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the turbulent energy is 



3M corc At£ 



81n(2) m 



(1) 



the mean turbulent Mach number is A4 — 0.79, and the 
ratio of turbulent to gravitational energy is 7 TURB = 0.010. 

We note that the core mass we are using is ~ 10 times 
higher than the value derived by Andre et al. (2007), who 
find M CORE = O.13M0 using 1.2 mm observations and as- 
suming a constant dust temperature of 20 K. On the other 
hand, it is ~ 3 times lower than the mass derived by |Sinu>| 
son et al.l d2008b, who find M c 



= 3.18M 



0- 



Finally, the 

virial mass for the given velocity dispersion as derived by 
Andre et all, ( |2007[ ) is M CORE = 0.71M Q . Given the signifi- 
ant spread in observational mass estimates for A-MM8 we 
have adopted the revised and intermediate value of Simp- 
son (20101. Our particular choice of parameters results in 
the core being highly supercritical, with a THERM + 7 T urb ~ 
0.027. Hence the results may only pertain to such highly su- 
percritical cores. We will explore the consequences of adopt- 
ing larger values of a THERM and/or 7 TURB , in a subsequent 
paper. 



method uses the density, p i , temperature, T i , and gravi- 
tational potential, tp t , of an SPH particle i to estimate a 
mean column-density, = Ti(p i ,tp i ), and a mean optical 
depth, f t = r(p i ,T i ,^p i ), through which the particle cools 
and heats. This optical depth includes contributions from 
dust, lines and free-free processes, and accounts for the vari- 
ation in the opacity in the cooler, less dense material that is 
presumed to surround particle i. The net radiative heating 
rate for the particle i is then 



dm 
dt 



4a SB {It -If) 



(2) 



The positive term on the right-hand side represents heat- 
ing by the background radiation field, and ensures that the 
gas and dust cannot cool radiatively below the background 
radiation temperature T Q , which we set to T Q = IK. The 
energy equation then takes account of compressional heat- 
ing, viscous heating, radiative heating by the background 
and radiative cooling. It has been extensively tested against 
detailed numerical ( Masunaga & Inutsuka 2000 ; Boss & Bo- 
denheimer||T979"l |Boss fc Myhill||1992l jWhitehouse fc Bate 



[20061 and analytical results (Spiegel 1957| |Hubeny||1990[ |, 
and performs well in both the optically thin and optically 
thick regimes. 



2.3 Gravity and hydrodynamics 



We use the SEREN SPH code ( |Hubber et aTpOU) , which 
is parallelized using OpenMP and designed for star forma- 
tion simulations. It has been extensively tested, and ap- 



plied to a wide range of problems (e.g. Stamatellos & Whit- 
worth|[20l0l |Stamatellos et al.||2l)TT| |Bisbas et al.||2009 
2010). It includes both the traditional SPH formulation 



| Monaghan||l992 1 and the more recent grad-h SPH formu- 
lation (Price & Monaghan 2004), which we use in this pa- 



per. To solve the SPH equations, we employ the symplectic 
2nd-order Leapfrog-KDK integrator, in conjunction with a 
block time-stepping scheme. We invoke additional features 
within the basic SPH algorithm, such as the Balsara vis- 



cosity switch (Balsara 19951 to reduce artificial shear vis- 



cosity. We use a Barnes-Hut octal-spatial tree (Barnes & 



Hut 1986) with the GADGET-style multipole acceptance 
criterion (Springel et al. 20011. Each SPH particle has a 



3 RESULTS 

3.1 Core angular momenta 

As fc MIN increases, the coherence lengths of the most ener- 
getic turbulent modes (< A MAX /2 ~ -R C ore/^min) decrease. 
Since these modes are uncorrelated, the specific angular mo- 
mentum, 



(3) 



which is compounded by contributions from all the different 
modes, also decreases in magnitude. Fig. [l] shows the vari- 
ation of j CORE with fc MIN at t = 0. For each value of fc MIN 
we simulate five different realizations by invoking five differ- 
ent seeds, and this produces an ~0.5dex spread in j CORE . 
However, there is an underlying systematic variation that, 



mass of m SPH = 10 M Q , resulting in a mass resolution of m the interval 1 < ^ < 4j can be approximated by 
10 M (a Jupiter mass). For random seed 500, we demon- 
strate convergence by re-simulating the fiducial setup using 
10 times as many SPH particles, i.e. with m SI 



x 10 19 cm 2 s" 



'k: 



(4) 



= 1O-°M . 

These runs are referred to as the '_hr' simulations. Gravi- 
tationally bound condensations that have a central density 



Pc > Psi 



10 gem are replaced with sinks, and sub- 



sequently grow by accretion, using a new algorithm (Hubber, 
Walch & Whitworth, 2011) that (a) ensures excellent numer- 
ical convergence, and (b) broadcasts the angular momentum 
of the accreted material to the surrounding gas (rather than 
assimilating it, which would be non-physical). 



2.4 Energy equation and radiative transfer 



At smaller fc MIN , j CORE approaches the maximum value con- 
sistent with the amount of turbulent energy (see below). In 
addition, much of the turbulent energy is invested in modes 
of such long wavelength that it results in bulk motion of the 
core, rather than intrinsic spin. The overall range of values 
is 10 17 - 8 cm 2 s" 1 < j CORE < lO 20 ' 2 cm 2 s" 1 . 

We note that for a critical Bonnor-Ebert sphere, the 
ratio of rotational to gravitational energy is given by /9 ROT — 
Ocore/^) 2 - where ^ = 0.644(GM ( 
MM8, j. = 1.6 x 10 21 cm 2 s" 1 and hence 



p \l/2 pi A 

COBE A COREi ■ rul A 



We use the radiative diffusion approximation of |Stamatellos| 
|et al.| ( |2007[ ) (RAD-WS method) to solve the energy equa- 
tion and evaluate radiative transfer effects. The RAD-WS 



(r 



J CORE 

x 10 21 cm 2 s" 



(5) 



Since /3 ROT ^ 7 T urb ~ 0-01 (^ ne rotational energy cannot 
exceed the turbulent energy), there is a maximum specific 
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Figure 1. Magnitude of the specific core angular momentum 
J'core as a function of fc MIN (lower abscissa), or equivalently 
^max/^core (upper abscissa). Results obtained with the same 
seed but different fc MIN are connected with dashed lines. The thick 
black line is the best fit to the data in the interval 1 < fc MIN < 4 
(see Eqn. B. 

angular momentum j MAX = 7itm B Ja ~ ^ x 10 20 cm 2 s _1 . 
Evidently this rather small j MAX is a consequence of small 
R CORE (observational estimates of j tend to derive from 
more extended cores) and small 7 TURB . 

3.2 Disk density distributions 

At t 50 , the time at which 50% of the initial core mass has 
been converted into protostars, we define a Cartesian frame 
of reference, (x IF ,y IF , z IF ), in which the z IF -axis is aligned 
with the largest principal moment of inertia of the remain- 
ing dense (p > 10 -12 gcm -3 ) material. The algorithm for 
doing this is described in Appendix A. Provided there is 
a single dominant primary disk, z IF is then aligned with 
its rotation axis. Fig. 2 displays a montage of false-color 
images of the density on the z IF — plane for the entire 
ensemble of simulations of A-MM8. Each column of images 
corresponds to a different value of k MIN , and each row to a 
different seed. From these plots we see that (a) the sizes of 
disks tend to decrease with increasing fc MIN ; and (b) that for 
fc MIN = 1 and 2, multiple protostars usually form, whereas 
for k MIN =1/2 and 4, only single stars are formed. 

3.3 Resolution study 

We re-simulate three setups for seed 500 (fc MIN = 
1/2, 1, and 4) with ten times higher resolution (run 500-hr), 
i.e. a total of 1,280,000 particles and m SPH = 10 -6 . The re- 
sults of these simulations are shown in Fig. [3] With regard 
to the density profiles and the disk radii, we find remarkably 
good convergence between low and high resolution runs (see 
also Table [oT2| . The disk masses are also in reasonable agree- 
ment, considering the fact that the end times are slightly 
different. For fc MIN = 1/2 we could not follow the simulation 
until 50% of the core mass collapsed into the sink because 
of CPU time limitations. Therefore, the runs are mistimed 
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Table 1. Properties of protostellar disks and protostars. Col- 
umn 1 gives fc MIN • Column 2 gives the seed used to generate the 
initial turbulent velocity field. The extension '_hr' denotes the 
high-resolution runs. Column 3 gives the time at which half the 
initial core mass has been converted into protostars, i 50 , in kyr. 
Column 4 gives the disk mass, M mSK , in M Q . Column 5 gives 
the disk radius, R DlSK , in AU. Column 6 gives the masses of the 
protostars formed, , in M Q . 



and the disk masses cannot be strictly compared in this case. 
For fc MIN = 1, the sink masses grow a bit quicker in the high 
resolution run (t 50 = 14.2 kyr instead of t 50 = 15.2 kyr) lead- 
ing to a smaller disk mass of O.14M0 rather than 0.31Mq. 
This difference may be caused by the different fragmentation 
properties as in run 500-hr only two instead of three sink 
particles form. Despite the different number of sink parti- 
cles, the disk density distributions and the 'system masses' 
within the two identifiable, individual condensations (see 
Fig. 2 and Fig. [3} are very similar, i.e. M t = 0.323 and 0.318 
for 500-hr and M t = 0.36 and (0.27 + 0.01) = 0.28 for 500. 
For fc MIN = 4, we find good agreement of all interesting 
quantities. Overall, the results of our study are only weakly 
dependent on resolution. 



3.4 Global disk properties and scaling relations 

There are many types of disk: circumstellar, circum-binary 
and circum-system (i.e. enclosing higher order systems). 
Here we focus our discussion on the primary accretion disk, 
which is the most massive and extended disk in the sim- 
ulation (at t 5Q ), and always surrounds the protostar with 
the highest final mass - although this is not always the first 
protostar to form. The typical primary disk is rather small, 
7? DISK ~ 10 to 30 AU, and has mass M DISK ~ 0.1 to 0.6 M Q . 
In general more extended disks are more massive, and the 
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Figure 2. Montage of false-color images of the density on the midplanes of the protostellar disks formed in the entire ensemble of 
simulations of A-MM8. Each column corresponds to a different value of fc MIN (from left to right: 1/2, 1, 2, and 4), and hence to a 
different A MAX = 2R CORB /fc MIN (from left to right: 4i? CORE , 2i? CORB , i? CORE , and iJ CORE /2). Each row corresponds to a different 
seed for generating the initial turbulent velocity field. The false-color encodes the same range of density on all plots: 10 — 15 gem -3 to 
10 -11 gem -3 . However, the linear sizes of the frames are different for different seeds, varying from 100 AU to 400 AU. Black dots mark 
the positions of sink particles. 
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Figure 3. Same as Fig. 2 but for the high resolution runs with seed 500 (500-hr). 




most extended disks display signs of gravitational instabil- 
ity. However, this only results in the formation of spiral arms 
and the transport of angular momentum; there are no pro- 
tostars formed by disk fragmentation. Table 1 lists for each 
simulation, the value of fe MIN , the seed used to generate the 
turbulent velocity field, the time to convert half the core 
mass into stars, t 50 , the mass and radius of the primary disk 
at this time, M DISK and -R DISK , and the masses of the pro- 



tostars formed. Fig. 3 shows the dependencies of -R DISK and 
M DISK on fc MIN and j CO REi an d Fig- 4 shows the correlation 
between M DISK and R DISK ■ The linear fits on these figures, 
and their uncertainties, obtained by \ 2 minimization, are 
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Figure 5. The correlation of disk mass, Af DISK , with disk radius, 
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4 DISCUSSION 

4.1 Filament formation 

Large-scale filaments play a critical role in core fragmenta- 
tion, and in the formation of the primary disk. First, large- 
scale filaments provide alternative sites (alternative to the 
center of mass of the core) where material can converge and 
form secondary protostars. Second, large-scale filaments de- 
liver large parcels of material with disparate angular mo- 
menta into the center, where they accumulate in a large 
disk around the primary protostar. Figs. [6] demonstrates 
that large-scale filaments are generated by turbulence with 
small fe MIN . 

The right-hand column of Fig. [6] shows false-color im- 
ages of the maximum density on lines of sight parallel to 
the z IF axis, at time t = 15kyr, for different values of k MTN . 
All the simulations presented in Fig. [6] derive from the same 
seed, but the results obtained with other seeds are statisti- 
cally similar. We see that for fe MIN = 1/2 the high-density 
gas (n > 10 6 H 2 cm~ 3 ) is concentrated in large-scale fila- 
ments. However, as fc MIN is increased, the strength and co- 
herence of the filaments declines, and by fc MIN = 4 there are 
no noticeable filaments. Filament formation can be under- 
stood in terms of the forces shaping the core. We neglect 
the pressure force, Fp, since Fp is initially smooth, and fo- 
cus on the ratio of centrifugal to gravitational force Fc/Fq. 
For each SPH particle, i, we compute F c = r; A Vi| 2 /|ri| 3 
and Fq — GAf(|r;|)/|r;| 2 , where and v< are the posi- 
tion and velocity of particle i relative to the centre of mass, 



Figure 7. Panel A: Mean (diamond) and variance (error bar) of 
the number of protostars formed, as a function of fcjvilN- Panel B: 
Mean (diamond), variance (error bar), and range (mandarine line) 
of the disk radii at t 50 as a function of the number of protostars 
formed. 



and M(|r;|) is the mass interior to radius |r,|. Parcels of 
gas with low centrifugal support, i.e. small Fc/Fq, collapse 
first, and neighboring parcels are then drawn into the re- 
gions they vacate, creating preferred accretion streams, i.e. 
filaments. The left-hand column of [6] shows false color im- 
ages of the minimum value of Fq/Fq found on each line of 
sight. For fc MIN =1/2, there are well defined structures with 
low Fc/Fq that can be related to the filaments illustrated 
in the corresponding right-hand image. However, as fc MIN is 
increased, structures with low Fc/Fq become increasingly 
small and incoherent. 

A complex of filamentary structures on scales of 1000 
AU, very similar to our case with small k Mm , has recently 
been observed in the envelopes of Class cores by |Tobin| 
|et al.| l |2010) using SPlTZER. [Tobin et al.| ( [2010[ ) note that this 
complex envelope structure is spatially distinct from possible 
outflow cavities, and explicitly suggest that it results from 
the collapse of prestellar cores with initial non-equilibrium 
structures. 



8 S. Walch, A. Whitworth, P. Girichidis 




Figure 6. Left COLUMN. The minimum value of the ratio of centrifugal to gravitational potential force, F^/Fq, found along the line of 
sight, for the initial conditions of the simulations with seed 200 and different A: MIN . For fc MIN = 1/2 we find coherent features, whereas 
the distribution of small Fq/Fq becomes increasingly random with increasing fe MIN . Right COLUMN. The maximum density along the 
line of sight after 15kyr of evolution for the simulations with seed 200 and different k M1N . Filaments have formed where Fq/Fq was 
initially small. 



The influence of the turbulent perturbation scale on prestellar core fragmentation and disk formation. 



4.2 Protostellar multiplicity 

From Table 1 it appears that cores with fc MIN = 1/2 and 4 
only spawn single stars, whereas cores with fc MIN = 1 and 2 
tend to spawn 2 or 3 stars. In other words, multiple systems 
are formed only if A MAX /ii OORE ~ 1 to 2. 

This is illustrated in Panel A of Fig. [7] where we plot 
the mean number of stars formed, and its variance, against 
k um . There are no instances of disk fragmentation. In all 
cases where secondary stars form they form by dynamical 
filament fragmentation. Even where two (or three) stars end 
up in a close binary (triple) system with a circum-binary 
(circum-system) disk, the components have always formed 
by filament fragmentation, with their own independent cir- 
cumstellar disks, and then fallen into the center and cap- 
tured one another. This result was already suggested by the 
results of Walch et al. (2010), but the current study places 
it on a quantitative footing. 

Panel B of Fig. 6, shows that by t 50 there is a wide 
range of disk sizes around single stars, but the disks around 
stars in multiple systems have been truncated by mutual 
tidal interactions. When k MIN = 1/2, the filaments do not 
fragment because they fall into the center very rapidly and 
consequently the material in them is stretched. On arrival in 
the center, much of this material is initially parked in a mas- 
sive extended disk around the primary protostar. However, 
despite being massive and extended, this disk does not frag- 
ment. Repeated perturbations due to irregular infall from 
the filaments (a) maintain a relatively high velocity disper- 
sion in the disk (and hence a high Toomre Q parameter), 
(b) excite density waves that transport angular momentum 
by gravitational torques, thereby facilitating accretion onto 
the central primary star, and (c) shear proto-condensations 
apart. In a similar vein, |Hayfield et al.| ( |2010l ) have recently 
shown that disks in binary systems tend to be more stable 
towards fragmentation than disks around single stars. 

When fc MIN = 4, there are no significant filaments, so 
dynamical fragmentation is suppressed. The disks that form 
around the primary protostar are too small to fragment. 

We stress that these results are for a specific low-mass, 
low-turbulence core. The critical value of kyim, above which 
fragmentation is inhibited, probably increases with increas- 
ing core mass, since a small low-contrast filament is more 
likely to fragment if it is more massive. Also, we might ex- 
pect more fragmentation for increased levels of turbulence. 



5 COMPARISON WITH PREVIOUS WORK 

There have been several other studies of the collapse and 
fragmentation of low-mass, turbulent cores. 

|Goodwin et al.| ([2004a 2004b) simulate the collapse and 
fragmentation of cores having mass M CORE =5.4 M Q and ra- 
dius R C ohe = 50, 000 AU, modeled with a Plummer-like den- 
sity profile. These are SPH simulations, using a barotropic 
equation of state and sink particles. Different levels of tur- 



bulence are considered, 0.01 ^ 7 T 



< 0.25, with P. ocfc" 



= 1 and purely solenoidal modes; for each case many re- 



alizations are performed. In contrast to the simulations pre- 
sented here (which result in fragmentation for 7turb 

= 0.01, 

provided 1<A; MIN <2), Goodwin et al. obtain fragmentation 
only when 7 TURB > 0.05. This is due to the fact that their 



cores have a much higher level of thermal support than ours, 
0.30<a THERM <0.45, so the gas is less readily compressed, 
and to the fact that solenoidal modes produce less compres- 
sion than the thermal mix of solenoidal and compressive 
modes that we use. 



In a related study Goodwin et al. (20061 extend the 
study with 7 TURB = 0.10 to different turbulent power spec- 
tra, P k <xk~ n with n = 3, 4, 5. They find that with higher 
n (i.e. a higher concentration of power at long wavelengths), 
there is more fragmentation, and the protostars formed have 
somewhat lower masses. 

These simulations have been repeated by Attwood et al. 



(20091, with the same initial conditions, but solving the en- 



ergy equation and treating the associated transport of cool- 
ing radiation, instead of using a barotropic equation of state. 
The main differences in the results are (i) that fragmentation 
is more efficient (larger numbers of protostars are formed), 
and (ii) that the binary systems have shorter periods, higher 
eccentricities, and smaller mass-ratios. 

|WalcheTaLl ( [20T0l ) use SPH to simulate the collapse and 
fragmentation of cores having mass Af CORE =6.1 M and ra- 
dius R CORE = 17, 000 AU, modeled as marginally supercrit- 
ical Bonnor-Ebert spheres (£ B = 6.9, density increased by 
10%). The cores have an initial ratio of thermal to gravita- 
tional energy a THERM =0.74, and are contained by an exter- 



nal pressure P E 



9 x 10" 11 erg cm -3 . 



The turbulent ve- 



locity field is characterized by a mean Mach number M = 1 
(i.e. transsonic turbulence) a power spectrum P k tx fc -4 , 
fc MIN = 1, and a thermal mix of solenoidal and compressive 
modes. A large ensemble of cores is generated, and from 
these a representative subset, having specific angular mo- 
menta spanning the range 0.1 < (j CORE /10 21 cm 2 s _1 ) < 2.7, 
is extracted and evolved. The energy equation is solved us- 



ing the molecular-line cooling rates of Neufeld et al. ( 1995 1. 
As in the simulations presented here, Walch et al. find that 
dynamical filament fragmentation dominates over disk frag- 
mentation. However, the gas in their simulations is much 
hotter (because dust cooling is not included), so the disks 
that form are more swollen. In addition, since they do not 
use sink particles, they are unable to follow the simulations 
to the point where multiple protostars with circumstellar, 
circum-binary and circum-system disks are formed. 



Offner et al. ( 2008 I simulate the collapse and fragmen- 



tation of turbulent cores, using AMR and a barotropic equa- 
tion of state. Their cores are produced in a large-scale sim- 
ulation of a collapsing molecular cloud, with either driven 
or decaying turbulence. Individual cores are then followed 
at higher resolution. | Offner et al.| ( |2008[ ) find that simula- 
tions with decaying turbulence form on average more low- 
mass protostars than simulations with driven turbulence. 
| Offner et al.| ([2009 ) simulate fragmentation in a turbulent 
box, using AMR and solving the energy equation. They show 
that radiative feedback from the forming protostars inhibits 
disk fragmentation, thereby reducing the number of low- 
mass multiple systems formed (see also Bate 2009c I. Further 
analysis of these results ( Offner et al. 2010 1 suggests that dy- 



namical filament fragmentation is the dominant mechanism 
forming low-mass stars and binary systems, rather than disk 
fragmentation. The material in filaments is sufficiently far 
from protostellar radiation sources to keep cool and frag- 
ment, whereas the material in disks around newly formed 
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protostars is close and gets heated up so that it does not 
fragment. 

However, 



Stamatellos et al. (20111 show that, if accre- 



tion onto a protostar is episodic (as is believed to be the 
case), the luminosity is also episodic, and the duty cycle has 
sufficiently long low-luminosity periods for the outer parts 
of a massive accretion disk to cool down and fragment. Disk 
fragmentation may therefore still be a viable mechanism for 
forming brown dwarfs. 



6 CONCLUSIONS 

We have performed an ensemble of SPH self-gravitating 
radiation-hydrodynamic simulations to demonstrate that - 
in low-mass turbulent cores - the largest wavelength in the 
turbulent spectrum has a critical bearing on the outcome 
of collapse and fragmentation. Specifically, if all other pa- 
rameters (the initial critical Bonnor-Ebert density profile, 



0.017, 7 T 



: 0.010, n = -d\n[P k ]/d\n[k] =4) 



are held fixed and A MAX is varied, 

• the mean specific angular momentum of the core in- 
creases approximately as the square of the the largest wave- 
length, j CORE ~ A^ AX , for l/2<A MAX /i? CORE <2; 

• filaments form in the regions where the centrifugal sup- 
port is weakest, and therefore the material collapses fastest; 

• the size and coherence of filaments therefore increases 
with increasing A MAX ; 

• dynamical filament fragmentation (cf. 



Offner et al. 



2010) is the dominant (only) fragmentation mechanism, 

• and hence fragmentation and multiple star formation 
only occur for 1/2<A MAX /R 

CORE r~> ^ ' 

• the primary (i.e. most massive) protostellar disk has 
mass and radius which scale approximately as M DISK ~ 

^MAX aI1C ^ ^"DISK ~~ \lAX 1 

• massive extended disks form (for large A MAX ) where the 
filamentary inflows deliver material with disparate specific 
angular momentum; 

• but these disks do not fragment, because the inflowing 
material maintains a large velocity dispersion and therefore 
the gravitational modes excited in the disk are only strong 
enough to redistribute angular momentum and facilitate ac- 
cretion onto the central protostar. 

The global parameters of a core do not completely spec- 
ify the initial conditions for a simulation. In particular, the 
initial turbulent velocity field is stochastic. Consequently 
there is considerable variance amongst different realizations 
of the same parameter set, and the conclusions listed above 
should be interpreted as statistical. 
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Radius [AU] 

Figure Al. Density profiles along the principal axes in the local 
inertial frame for the disk formed in the simulation with fc MIN = 
1/2 and seed 400. The dotted vertical lines mark the radii at 
which the density along the x IF - and j/ IF -axes drops below the 



threshold density of p, 



THRESH 



10~ 



APPENDIX A: DISK DEFINITION 

In order objectively to identify disks and e.auxvaluate their 
masses and radii, we apply the following procedure. First, 
we isolate all the material with density p > 10~ 12 gcm -3 . 
Second, we compute the moment of inertia tensor for this 
material, and thereby define a new local co-ordinate system 
using the principal axes of inertia, x IF , y lF , z IF (where the 
subscript "if" stands for inertial frame). The properties of 
the dense material are now analyzed relative to this new co- 
ordinate system. In particular, z IF is allocated to the largest 
principal moment of inertia, and therefore, if the material is 
in a disk, z 1F is its rotation axis. 

In order to ascertain whether there is a disk, we com- 
pute the logarithmic density profile along each of the axes 
of inertia, and smooth these profiles using a box-car aver- 
aging technique. For a disk, the profiles along a; IF and y 1F 
are very similar to one another, and the third, along z IF is 
significantly steeper and less extended. 

Because disks are embedded in, and grow from.aux, fil- 
aments, we locate the edge of the disk at the first point 
along the a; IF - and j/ IF -axes where the density is below 



3 and the second derivative of the 

2i „/j/l„„ „\2 



Pthresh = 1U gem 

logarithmic density profile is zero, d 2 log 10 p/d(log lQ r) z = 
Fig. Al shows the profiles along the principal axes for the 
disk formed in simulation with fc MIN = 1/2 and seed 400. 
This is the disk illustrated in the left column of the middle 
row of Fig. 2. All the disks in Fig. 2 have been identified 
in this way, and are viewed face-on down z IF in their local 
inertial frame. 

The dotted vertical lines on Fig. Al mark the radii at 
which the density first falls below P THR esh > at §8 AU on x IF , 
and at 95 AU on y 1F , respectively. We identify the edge of the 
disk where the second derivative of the logarithmic density 
next falls to zero. This gives a mean radius of -R DISK = 94 AU, 
in this case. The density typically drops very steeply inside 
this radius, and therefore the resulting estimate of the disk 
mass (M DISK =0.57M„ in this case) is robust. 
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